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1 Introduction 



Standard derivative pricing theory is based on the assumption of frictionless 
and perfectly liquid markets. In particular, it is assumed that all investors are 
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small relative to the market so that they can buy arbitrarily large quantities 
of the underlying assets without affecting its price (perfectly liquid markets). 
Given the scale of hedging activities on many financial markets this is clearly 
unrealistic. Hence in recent years a number of approaches for dealing with 
market illiquidity in the pricing and the hedging of derivative securities have 
been developed. The financial framework that is being used differs substan- 
tially between model classes. However, as shown in Section 2 below, in all 
models derivative prices can be characterized by fully nonlinear versions of 
the standard parabolic Black-Scholes PDE. Moreover, these nonhnear Black- 
Scholes equations have a very similar structure. This makes these nonlinear 
Black-Scholes equations a useful reference point for studying derivative asset 
analysis under market illiquidity. In Sections 3 we therefore study scaling prop- 
erties and compute explicit solutions for a typical fully nonlinear Black-Scholes 
equation; Section 4 is devoted to studying properties of solutions and sensi- 
tivities with respect to model parameters. The relevant literature is discussed 
in the body of the paper. 



2 Illiquid mcirkets and nonlineeir Black-Scholes equations 

In order to motivate the subsequent analysis we present a brief synopsis of 
three different frameworks for modeling illiquid markets. We group them under 
the labels a)transaction-cost models] b) reduced-form SDE-models; c) reaction- 
function or equilibrium models. In particular, we show that the value function 
of a certain type of selffinancing strategies (so called Markovian strategies) 
must be a solution of a fully nonlinear version of the standard Black-Scholes 
equation. In all models there will be two assets, a risk-free money-market 
account B which is perfectly liquid and a risky and illiquid asset S (the stock) . 
Without loss of generality, we use the money market account as numeraire; 
hence Bt = 1, and interest rates can be taken equal to zero. 

(Quadratic) transaction-cost models. The predominant model in this 
class has been put forward by Cetin, Jarrow and Protter [7]. In this model 
there is a fundamental stock price process following geometric Brownian 
motion. 



for constants G M, o" > and a standard Brownian motion W. The trans- 
action price to be paid at time t for trading a shares is 



dS'^ = i^S^dt + aS^dWt 



(2.1) 




a eR, p> 0, 
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where p is a liquidity parameter. 

Intuitively, in the model (2.2) a trader has to pay a spread whose size rela- 
tive to the fundamental price, Sf{e''°' — 1), depends on the amount a which is 
traded. As shown by [7] , this leads to transaction costs which are proportional 
to the quadratic variation of the stock trading strategy. In order to formu- 
late this statement more precisely we consider a selffinancing trading strategy 
{^t,Vt)t>o giving the number of stocks and the position in the money market 
for predictable stochastic processes $ and rj. The value of this strategy at time 
t equals Vt = ^tSt + Vt- 

Consider now a Markovian strategy, i.e. a trading strategy of the form $t = 
S^) for a smooth function 0. In this case $ is a semimartingale with 
quadratic variation given by 

t 

= j (<t>s{s,Sl)aSf)'ds- 



see for instance Chapter 2 of [22] . Then Theorem A3 of [7] yields the following 
dynamics of Vt 

dVt - S^)dS^ - pS^d^t = S^)dS^ - pS^ (Mt, S^)aS^y dt.{2.3) 

Note that according to (2.3), the market impact of the large trader leads to 
an extra transaction-cost term of size p /q d[^]s in the wealth dynamics. 

Suppose now that u and are smooth functions such that u{t, S^) gives the 
value of a selffinancing trading strategy with stock position 0(t, S^). According 
to the Ito formula, the process {u{t, S^))t>o has dynamics 

du{t, S^) = us{t, S^)dS^ + (ut{t, S^) + ^a'{S^fuss{t, 5°)) dt . 

Comparing this with (2.3), using the uniqueness of semimartingale decomposi- 
tions it is immediate that u must satisfy the PDE ut+^a'^S^uss+pS^a^^'s — 
and that moreover = us- The last relation gives (ps — uss, so that we obtain 
the following nonlinear PDE for u 

ut + ^a^S'uss (1 + 2pSuss) = . (2.4) 

When pricing a derivative security with maturity date T and payoff /i(S't) 
for some function h : [0, oo) — > M we have to add the terminal condition 
u{T,S) = h{S), S > 0. For instance, in case of a European call option with 
strike price K > we have h{S) = maxj^" — K, 0}. 
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Note that the original paper [7] goes further in the analysis of quadratic trans- 
action cost models. To begin with, a general framework is proposed that con- 
tains (2.2) as special (but typical) case. Moreover, conditions for absence of 
arbitrage for a general class of trading strategies - containing Markovian trad- 
ing strategies as special case - are given, and a notion of approximative market 
completeness is studied. 

Reduced-form SDE models. Under this modeling approach it is assumed 
that investors are large traders in the sense that their trading activity affects 
equilibrium stock prices. More precisely, given a liquidity parameter p > and 
a semimartingale $ representing the stock trading strategy of a given trader, 
it is assumed that the stock price satisfies the SDE 

dSt = aStdWt + pStd^t ■ (2.5) 



The intuitive interpretation is as follows given that the investor buys (sells) 
stock (A$4 > 0) the stock price is pushed (downward) upward by pSt-Z^^t] 
the strength of this price impact depends on the parameter p. Note that for 
p = the asset price simply follows a Black- Scholes model with reference 
volatility cr. The model (2.5) is studied among others in [10], [11], [14] or [18]. 

In the sequel we denote the asset price process which results if a large uses 
a particular trading strategy $ by 5*. Suppose as before that the trading 
strategy is Markovian, i.e. of the form = 0(t, St) for a smooth function 
0. Applying the Ito formula to (2.5) shows that S'^ is an Ito process with 
dynamics 

dSf = v'f'it, St)StdWt + h^{t, Sf)Stdt 
with adjusted volatility given by 

1 - pS(ps{t,S) 



see for instance [10] for details. In the model (2.5), a portfolio with stock 
trading strategy $ and value V is termed selffinancing, if satisfies the equation 
dVt — ^tdSf. Note that the form of the strategy $ affects the dynamics of S; 
this feedback effect will give rise to nonlinearities in the wealth dynamics as we 
now show. Suppose that Vt = u{t, St) and ^t = 4>{t, St) for smooth functions u 
and 0. As before, applying the Ito formula to the process {u(t, Sf))t>o yields 
(pt = Us- Moreover, u must satisfy the relation ut + ^{v'^)'^{t, S)S'^uss — 0. 
Using (2.6) and the relation (ps — uss we thus obtain the following fully 
nonlinear PDE for u{t, S) 
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Again, for pricing derivative securities a terminal condition corresponding to 
the particular payoff at hand needs to be added. 

Equilibrium or reaction-function models. Here the model primitive is 

a smooth reaction function ip that gives the equilibrium stock price St at 
time t as function of some fundamental value Ft and the stock position of a 
large trader. A reaction function can be seen as reduced-form representation 
of an economic equilibrium model, such as the models proposed in [12], [21] 
or [23]. In these models there are two types of traders in the market: ordinary 
investors and a large investor. The overall supply of the stock is normalized 
to one. The normalized stock demand of the ordinary investors at time t is 
modelled as a function D{Ft,x) where x is the proposed price of the stock. 
The normalized stock demand of the large investor is written the form p^fi 
p > is a parameter that measures the size of the trader's position relative 
to the total supply of the stock. The equilibrium price St is then determined 
by the market clearing condition 

D{Ft,St)+p<^t = l. (2.8) 



Under suitable assumptions on D equation (2.8) admits a unique solution, 
hence St can be expressed as a function ip of Ft and i.e. St = ip{Ft, p^t)- 
For instance we have in [21] '^{f,a) — f exp{a); the model used in [12] and 
[25] leads to the reaction function ip{f, a) — f/{l — a) . The reaction-function 
approach is also used in [15] and in [9]. 

Now we turn to the characterization of selffinancing hedging strategies in 
reaction-function models. Throughout we assume that the fundamental-value 
process F follows a geometric Brownian motion with volatility cr (see (2.1)). 
Moreover, we assume that the reaction function is of the form ip{f, a) = fg{a) 
for some increasing function g. This holds for the specific examples introduced 
above and, more generally, for any model where D{f, x) — U (f/x) for a strictly 
increasing function U : (0, oo) — > IR with suitable range. 

Assuming as before that the normalized trading strategy of the large trader is 
of the form p0(t, S) for a smooth function we get from Ito's formula, since 
St^g{p<f>{t,St))Ft, 

dSt = g{p(p{t, St)) dFt + pFtgaip^t, St))Mt, St) dSt + 6(t, St) dt; 

the precise form of b(t, St) is irrelevant for our purposes. Rearranging and 
integrating the expression (1 — pFtga{p(t){t, St))(f)s{t, St))~^ over both sides of 
this equation gives the following dynamics of S: 

dSt = , (rStdWt + b(t, St)dt . (2.9) 

^-p'MMstMt,St) 
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Remcirk 2.1 A very general analysis of the dynamics of seljfinancing strate- 
gies in reaction-function models (and generalizations thereof) using the Ito- 
Wentzell formula can be found in [Ij. 



A similar reasoning as in the case of the reduced-form SDE models now gives 
the following PDE for the value function u{t,S) of a selfhnancing strategy 

1 cr^ 

^*t+o7 -^S^ss^O. (2.10) 



In particular, for g{a) — exp{a) we have g — Qa and (2.10) reduces to equation 
(2.7); for g{a) = 1/(1 - a) as in [12], [25], we get the PDE 



Ut 



+ ^ 2^ '^SS = . (2.11) 

2(1- pus - pSuss) 



Nonlinear Black-Scholes equations. The nonlinear PDEs (2.4), (2.7) , 
(2.10) and (2.11) are all of the form 

Ut + ^(^^v{S, pus, puss)S^uss = , (2.12) 



where v{S,0,0) — 1. Since p is often considered to be small, it is of interest 

to replace v with its first order Taylor approximation around p = 0. It is 
immediately seen that for the equations (2.7) and (2.11) this linearization is 
given by v{S, pus, pUugs) ~ 1 + pSuss] replacing v with this first order Taylor 
approximation in (2.7) and (2.11) thus immediately leads to the PDE (2.4). 

Note that (2.12) is a fully nonlinear equation in the sense that the coefficient 
of the highest derivative is a nonlinear function of this derivative. A similar 
feature can be observed for the limiting price in certain transaction cost models 
under a proper rescaling of transaction cost and trading frequency; see for 
instance [2] . Nonlinear PDEs for incomplete markets obtained via exponential 
utility indifference hedging such as [3] or [19] on the other hand are quasi-linear 
equations in the sense that the highest derivative uss enters the equation in 
a linear way, similar to the well-known reaction-diffusion equations arising 
in physics or chemistry. Prom an analytical point of view the nonlinearities 
arising in both cases are quite different. In the latter quasi-linear case we have 
to do with a regular perturbation of the classical Black-Scholes (BS) equation 
but in the case (2.10) we have to do with singular perturbation, meaning 
that the highest derivative is included in the perturbation [16]. In case of a 
regular perturbation it is typical to look for a representation of a solutions to 
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a quasi-linear equation in the form 

TV 

u {S, t) = UBs{S, t) + Y. puliS, t) + 0(p^+^). (2.13) 

n=l 

For many forms of nonlinearities the uniform convergence as p — of the 
expansions of type (2.13) can be estabhshed. If the perturbation is of the 
singular type, the asymptotic expansion of type (2.13) typically breaks down 
for some S,t and some > 0. It can as well happens that the definition 
domain for the solutions of type (2.13) include just one point or that it is 
empty. It is very important to obtain explicit solutions for such models because 
in these case we can not hope to get good approximations for solutions by 
expansions of the type (2.13) in the whole region. 



3 Invariant Solutions for a Nonlinear Black-Scholes Equation 



Our goal in this section is to investigate the nonlinear Black-Scholes equation 
(2.7) using analytical methods. As we have just seen, this equation is typi- 
cal for the nonlinearities arising in pricing equations for derivatives in illiquid 
markets in generalQ Using the symmetry group and its invariants the partial 
differential equation (2.7) can be reduced in special cases to ordinary differ- 
ential equations. In [6] a special family of invariant solutions to the equation 
(2.7) was studied; in particular, the explicit solutions were used as test case 
for various numerical methods. In the present paper we study all invariant 
solutions to equation (2.7). 



3.1 General Results 



In order to describe the symmetry group of (2.7) we use results from Bordag 
[5]. In that paper the slightly more general equation 



with a continuous function A : (0, oo) (0, oo) is studied (note that for A = 1 
this equation reduces to (2.7)) and the following theorem was proved. 



^ Our approach can be applied to equations (2.4) and (2.11) as well; see [4,24] for 
details. 
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Theorem 3.1 The action of the symmetry group Ga of (3.14) ^^^^ o.iT' arbi- 
trary function X{S) is given by 



S = S, (3.15) 

i=t-\-a2e, (3.16) 

u — u-\-asSe-\-a4e, ee(— oo,cxd). (3-17) 

// the function X{S) has the special form X{S) — ujS^, k,u! e M, equa- 
tion (3.14) takes the special form 

2 c2 



In that case the symmetry group Ga has the richer structure 



S = Se'''% ee (-00,00), (3.19) 
t — t + a2e, 

+ - 1)' M 0, M 1 (3.20) 

ai(l — /ej 

u^ue"'' + asSee"'' + —(e"'' - 1), k^ 0, 

^2 = + ^5(e"i^ - 1) + 046, A; = 1. (3.21) 
ai 



The proof is based on methods of Lie point symmetries, i.e. the Lie symmetry 
algebras and groups to the corresponding equations were found; see [5] for 

details and [17,20,13] for a general introduction to the methodology. As we 
will see later, the solutions found in [5] for the case \{S) = S can be used to 
obtain the complete set of invariant solutions to equation (2.7). 

For X{S) — LoS'^ with A; = 0, 1 we can use (3.19) to obtain a new independent 
invariant variable z and (3.21) to obtain a new dependent variable v; these 
variables are given by 



z — log S -\- at, a ^ 0, 

v^uS^''-^\ A; = 0,1. (3.22) 
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Using these variables we can reduce the partial differential equation (3.18) to 
the ordinary differential equation 

'"^ + ^Ti n ^^2 = 0' 3.23) 

(1 - b{v^^ + ^v^)y 



where 9 = f^, a ^ 0,b = up ^ 0,^ = (— 1)*^, k = 0, 1. It is a straightforward 
consequence of Theorem 3.1 that these are the only nontrivial invariant vari- 
ables. In the sequel wc will determine explicit solutions for (3.23) for the case 
A; = which corresponds to the original equation (2.7). 

Remark 3.1 Equation (3.23) with k — and k — 1 are related to each other. 

Note however, that the relation between the corresponding solutions is not so 
straightforward, because (3.23) is nonlinear and we need real valued solutions. 
Hence the results from [5] - where (3.23) was studied for the case k = 1 - do 
not carry over directly to the present case k = 0, so that a detailed analysis of 
the case k — is necessary. 

To find families of invariant solution we introduce a new dependent variable 
y{z)^v,{z). (3.24) 



If we assume that the denominator of the equation (3.23) is different from 
zero, we can multiply both terms of equation (3.23) by the denominator of the 
second term and obtain 

yyl + 2e {y' - C + yz + {y' - + ) y = o, M 0.(3.25) 



We denote the left hand side of this equation by F{y, y^). The equation (3.25) 
can possess exceptional solutions which are the solutions of a system 

^^^1^ = 0, F(y,y,)=0. (3.26) 

oyz 



The first equation in this system defines a discriminant curve which has the 
form 

y(^) = I- (3.27) 

If this curve is also a solution of the original equation (3.25) then we obtain 
an exceptional solution to equation (3.25). We obtain an exceptional solution 
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if g = —4^, i.e. a — ^cr'^/S. It has the form 



(3.28) 



This solution belongs to the family of solutions (3.30) by the specified value 
of the parameter q. In all other cases the equation (3.25) does not possess any 
exceptional solutions. 

Hence the set of solutions of equation (3.25) is the union of the solution sets 
of following equations 




y = 0, (3.29) 
y={^±^)/b, (3.30) 

y^-[-Cy' + ly-:^,-\lrAiz-y]]-.y^^ (3-3i) 
y^ - [-^y' + ly-7h + \lrAij:-y]]7.y^^ (3-32) 

where one of the solutions (3.30) is an exceptional solution (3.28) by g = —4,^ 
for A: = 1. In case A; = we have ,^ = 1 so that solutions of (3.30) arc complex 
valued functions. We denote the right hand side of equations (3.31), (3.32) by 
f{y). The Lipschitz condition for equations of the type yz = f{y) is satisfied 
in all points where the derivative ^ exists and is bounded. It is easy to see 
that this condition will not be satisfied by 

y = 0, y = ^> y^oo. (3.33) 



Hence on the lines (3.33) the uniqueness of solutions of equations (3.31), (3.32) 
can be lost. We will study in detail the behavior of solutions in the neighbor- 
hood of the lines (3.33). For this purpose we look at the equation (3.25) from 
another point of view. If we assume now that z, y, yz are complex variables 
and denote 

y(^)=C, yz{z)^w, C,'weC, (3.34) 



then the equation (3.25) takes the form 

F(C, w) = C-^ + 2e (C^ - e + n, + (C^ - + ^) C = 0, (3 

where b ^ 0. The equation (3.35) is an algebraic relation in and defines 
a plane curve in this space. The polynomial F{(,w) is an irreducible polyno- 
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mial if at all roots Wr{z) of F{(^,Wr) either the partial derivative F(^{(,Wr) or 
Fy^{(, Wr) are non equal to zero. It is easy to prove that the polynomial (3.35) 
is irreducible. 

We can treat equation (3.35) as an algebraic relation which defines a Rieman- 
nian surface F : F{(,w) =0 oi w = w{() as a compact manifold over the 
(^-sphere. The function w{() is uniquely analytically extended over the Rie- 
mann surface F of two sheets over the sphere. We find all singular or branch 
points of w{() if we study the roots of the first coefficient of the polynomial 
F{(,w), the common roots of equations 

F{C,w)=0, F^{C,w) = 0, CweCUoo. (3.36) 



and the point ^ = oo. The set of singular or branch points consists of the 
points 

Ci = 0, C2 = ^, C3 = oo. (3.37) 



As expected we got the same set of points as in real case (3.33) by the study 
of the Lipschitz condition but now the behavior of solutions at the points is 
more visible. 



The points (2, (3 are the branch points at which two sheets of F are glued on. 
We remark that 

-(G)^-i(l-4)-*^ + -. (3.38, 



where t is a local parameter in the neighborhood of (2- For the special value 
oi q = 4^, i.e. k=0, the value w{(2) is equal to zero. 

At the point ^3 = 00 we have 

where t is a local parameter in the neighborhood of ^3. At the point Ci = 
the function w{Q has the following behavior 



q 1 

w{() ~ C ~^ Ci = 0, on the principal sheet, (3.39) 

"^^(0 K + ^ j C, C ^ Ci = 0, g 7^ 1, on the second sheet, (3.40) 
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w{C) ~ -2fe^C^ C ^ Ci = 0, g = 1, on the second sheet. (3.41) 



Any solution w{C) of an irreducible algebraic equation (3.35) is meromorphic 
on this compact Riemann surface F of the genus and has a pole of the order 
one correspondingly (3.39) over the point Ci = and the pole of the second 
order over ^3 = oo. It means also that the meromorphic function w{Q cannot 
be defined on a manifold of less than 2 sheets over the ( sphere. 

To solve differential equations (3.31) and (3.32) from this point of view it is 
equivalent to integrate on F a differential of the type and then to solve 

an Abel's inverse problem of degenerated type 



/ 



=z + const. (3.42) 



^(0 



The integration can be done very easily because we can introduce a uniformiz- 
ing parameter on the Riemann surface F and represent the integral (3.42) in 
terms of rational functions merged possibly with logarithmic terms. 

To realize this program we introduce a new variable (our uniformizing param- 
eter p) in the way 



46(p+l) ^ ^ 

Then the equations (3.31) and (3.32) will take the form 



f^f+^y^f + const, (3.45) 



(p-iMv + if + V;) 

p{p — l)dp 

(p+iMp-ir+^o 



— z + const. (3.46) 



The integration procedure of equations (3.45), (3.46) gives rise to relations 
defining a complete set of first order differential equations. In order to see 
that these are first order ordinary differential equations recall that from the 
substitutions (3.34) and (3.24) we have 



p^^l-%,. (3.47) 



We summarize all these results in the following theorem. 
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Theorem 3.2 The equation (3.23) for arbitrary values of the parameters 
q,b can be reduced to the set of first order differential equations which 
consists of the equations 

v,^0, v,^{l±2)/b, (3.48) 



and equations (3.31), (3.32) where y is defined by the substitution (3.24). The 
complete set of solutions of the equation (3.23) coincides with the union of 
solutions of these equations. 

To solve equations (3.45), (3.46) exactly we have first to integrate and then 
invert these formulas in order to obtain an exact representation of p as a func- 
tion of z. If an exact formula for the function p = p{z) is found we can use the 
substitution (3.47) to obtain an expUcit ordinary differential equation of the 
type Vz{z) = f{z) or another suitable type; in that case it is possible to solve 
the problem completely. However, for an arbitrary value of the parameter q 
inversion is impossible, and we have just an implicit representation for the so- 
lutions of the equation (3.23) as solutions of the imphcit first order differential 
equations. 

3.2 Exact invariant solutions for a fixed relation between S and t 

For a special value of the parameter q, namely for q = —4, we can integrate 
and invert the equations (3.45) and (3.46). For q = —4 the relation between 
the variables S and t is fixed in the form 

2 

z ^ log S-^t, (3.49) 

equation (3.45) takes the form 

(p-l)^(p + 2) = 2c exp(-3z/2), (3.50) 

and correspondingly the equation (3.46) becomes 

(p + l)2(p-2) = 2 c exp(-3^/2), (3.51) 



where c an arbitrary constant. It is easy to see that the equations (3.50) and 
(3.51) are connected by a transformation 

p —p, c — > — c. (3.52) 
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This symmetry arises from the symmetry of the underhning Riemann surface 
r (3.35) and corresponds to a change of the sheets on F. Using these symmetry 
properties we can prove following theorem. 

Theorem 3.3 The second order differential equation 

- "^TT^ir^ir^ = 0' (3.53) 

(1 - b{v^^ + v^)y 



is exactly integrable for any value of the parameter b. The complete set of 
solutions for b ^ is given by the union of solutions (3.57), (3.59) -(3.62) 
and solutions 

3 1 

v{z) — v{z) — -z -\- d, v{z) — —-Z + d, (3.54) 



where d is an arbitrary constant. The last solution in (3.54) corresponds to 
the exceptional solution of equation (3.25). 

For 6 = equation (3.53) is linear and its solutions are given by v{z) ~ 
di + d2exp — (32;/4), where di, d2 are arbitrary constants. 

Proof. Because of the symmetry (3.52) it is sufficient to study either the 
equations (3.50) or (3.51) for c G M or both these equations for c > 0. The 
value c = can be excluded because it complies with the constant value of 
p{z) and correspondingly constant value of Vz{z), but all such cases are studied 
before and the solutions are given by (3.54). 

We will study equation (3.51) in case c e M \ {0} and obtain on this way the 
complete class of exact solutions for equations (3.50)-(3.51) and on this way 
for the equation (3.53). 

Equation (3.51) for c > has a one real root only. It leads to an ordinary 
differential equation of the form 



Vz{z) = - 1 + 1 + ce-- + V2ce 



2 + e ^ ^ 



+ (l + ce ^■2+y2ce ^ + c^e-^'' 



, c> 0. 



(3.55) 
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Equation (3.55) can be exactly integrated if we use an Euler substitution and 
introduce a new independent variable 



r = 2(l + ce^ + V2ce^ +c2e-3M . (3.56) 



The corresponding solution is given by 



_ 2 

Vr{z) = -- (l + ce-^ + \/2ce-^ +c2e-3^) ' (3.57) 

^ 2 

log I ^1 + ce"^ + Vsce"^ + e-^^^ 



where d e R is an arbitrary constant. If in the right hand side of equation 
(3.51) the parameter c satisfies the inequality c < and the variable z chosen 
in the region 

z e f-^ln^,ooj (3.58) 



then the equation on p possesses maximal three real roots. These three roots 
of cubic equation (3.51) give rise to three differential equations of the type 
Vz = —(1 — p^{z))/b. The equations can be exactly solved and we find corre- 
spondingly three solutions Vi{z), i = 1, 2, 3. 

The first solution is given by the expression 



vAz) = — - — - cos arccos fl — IcI e 2^)^ (3.59) 
V3 ^ w 

4 , . fl 



— — log ^1 + 2 cos ^- arccos {l — |c| e 2 



— ^ log fsin (- arccos (^1 — |c| e 2 j'j'j + (i, 
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where d e R is an arbitrary constant. The second solution is given by the 
formula 

V2\z) — — - — - COS ( -TT — - arccos \ \ — |c| e 2 j 1 (3.60) 
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4 / /II / , , ^ 

— log —1 + 2 cos -TT arccos 1 — c e 2 

36 \ \3 3 ^ 

^ log (sin Qtt - - arccos (l - |c| e^))) + d, 



where 0? G M is an arbitrary constant. The first and second solutions are defined 
up to the point z = — | In |^ where they coincide (see Fig. 1). 

The third solution for 2; > — | In ||| is given by the formula 



, z 2 /2 2 / I I -3z \ \ 

V3,2{z) ^ (3^ + 3 arccos [l-\c\e 2 j 1 (3.61) 

log —1 + 2 cos -TT H — arccos ( — 1 + c e 2 

.36 V V .3 .3 V y 



36 ''V V3 
- log (cos (^-TT + - arccos (-1 + |c| e^))) + d, 

where G M is an arbitrary constant. In case 2; < | In j|| the polynomial (3.51) 
has a one real root and the corresponding solution can be represented by the 
formula 



,i{z) — —^ — ^ cosh ^^arccosh 1 + |c| e^)^ (3.62) 
— ^ log ^cosh ^-arccosh (—1 + |c| ^^~)^^ 

log ^—1 + 2 cosh ^^arccosh (—1 + |c| e^^^^ + d. 



4 

36 



The third solution is represented by formulas ^3,2 (2;) and ^3,1 (2;) for different 
values of the variable z. 

One of the sets of solutions (3.57), (3.59) -(3.62) for fixed parameters 6, c, d is 
represented in Fig. 1. The first solution (3.55) and the third solution given by 
both (3.61) and (3.62) are defined for any values of z. The solutions vi{z) and 
V2{z) cannot be continued to the left after the point z — —^\n ^ where they 
coincide. 

2 

If we keep in mind that z — log 5" — and u{S, t) = S v{z) we can represent 
the complete set of exact invariant solution of equation (2.7). The solution 
(3.57) gives rise to an invariant solution Ur{S,t) in the form 



Ur{S, t) = -^S (1 + cS-h^' + V2c5-fe^* + c2 S-^e^'' 
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Fig. 1. Plot of the solution v{z) given in (3.57) (thick solid line); vi{z) 
from (3.59) (short dashed line); V2{z) from (3.60) (long dashed line) and 
the third solution t'3,i(-z), i'3,2(-2) from (3.61), (3.62) (thin solid line). The 
parameters takes the values |c| = 0.35, q = —4, d = 0,b = 1. and the 
variable z G (—5,4.5). 



2 



-J_sll + cS-h"^' + \l2cS-ie^' + S'^e'tA ' (3.63) 



uop 



1 



+ (^l + c5-^eTr* + V2c5-ie%* + c2 5-3e"l-*J - 2 j + 5 + cis 

where d, (i2 G c > 0. This solution was obtained and studied in [6]. We 
describe now other invariant solutions from the complete set of invariant so- 
lutions. 

In case c < we can obtain correspondingly three real solutions if 



The first solution is represented by 



uxiS.t) = — -S I log^ - — 1 1 - — 5cos ( ^ arccos ( 1 - IcI S~^e^^ 
ijjp \ Q J up V3 ^ 
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Fig. 2. Plot of solutions Ur{S,t),ui{S,t),U2{S,t) and U3^i{S,t),U3^2{S,t) 
for the parameters a = 0.4, |c| = 0.5, q = —4,b = 1.0, d = 0. The 
variables S,t lie in intervals S S (0.,5.) and t £ [0,0.5]. All invariant 
solutions change slowly in t-direction. 



^ 5 log (^1 + 2 cos arccos (^1 - |c| S'te^*^ ^ ^ (3.65) 
Slog (sin (- arccos (l — \c\ S'^e^''^^^ + d S + d2, 



where d, (^2 € M, c < 0. The second solution is given by the formula 



U2(S,t) = S f logS* 1] S'cos ( -7T H — arccos ( —1 + IcI S ae^is * 

top \ 8 J tup V3 3 



^ 5 log f 1 + 2 cos (In + I arccos ( -1 + |c| S'h^^ ) ) ) (3.66) 



3(jjp \ \3 3 

S'log fsin (-TC + - arccos ( -1 + IcI S~'h^^ ]]] + d S + d2. 



Stop V V6 6 

where d, c?2 G M, c < 0. The first and second solutions are defined for the 
variables under conditions (3.64). They coincide along the curve 
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and cannot be continued further. The third solution is defined by 



us,2{S,t) = -^^S (^log S-^t^ (3.67) 
b COS - arccos — 1 + c Lb 2e le 



cup v3 
4 „, / n 



Scup 
16 



5'log ^-1 + 2 cos ^- arccos ^-1 + |c| 5" ae^"* * 
5* log ^cos Q arccos ^-1 + |c| ^'^e^*^^^ + d S + d2, 



Sup 

where d, d2 ^'M. and S,t satisfied the condition (3.64). 

2 

In case < 5" < (^^^ ^ exp (^^) the third solution can be represented by the 
formula 



U3,i{S, t) = -—S (logs - ^t] (3.68) 

LUp \ tS J 



2 /2 / 3 3<t2 

— 5" cosh I -arccosh I — 1 + IcI S~^e~* 
cup V3 V 

y^log fcosh f -arccosh f — 1 + IcI 5"" 
Sup V V6 V 



2e 16 



4 



S'log ^-1 + 2 cosh Qarccosh (^-1 + \c\ S 



3 Lop 
+ d S + d2, 

where d, d2 G R. The solution u{S,t) (3.63) and the third solution given by 
1*3,1) '^3,2 (3. 67), (3. 68) are defined for all values of variables t and 5" > 0. They 
have a common intersection curve of the type S — const, cxp(^t). The typical 
behavior of all these invariant solutions is represented in Fig. 2. 

The previous results we summarize in the following theorem which describes 
the set of invariant solutions of equation (2.7). 

Theorem 3.4 The invariant solutions of equation (2.7) can be defined by the 
set of first order ordinary differential equations listed in Theorem 3.2. 
If moreover the parameter q — —A, or equivalent in the substitution (3.22) 
we chose a = — o"^/8 then the complete set of invariant solutions of (2.7) can 
be found exactly. This set of invariant solutions is given by formulas (3.63)- 
(3.68) and by solutions 



u{S,t)^dS, u{S,t)^jS{logS-a^l), 

b 8 
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where d is an arbitrary constant. This set of invariant solutions is unique up 
to the transformations of the symmetry group Ga given by Theorem 3.1. 

The solutions Ur{S,t), (3.63), Ui{S,t), (3.65), U2{S,t), (3.66), 1x3 i(S', t) ,(3.67), 
'<J'3,2{S,t) (3.68), have no one counterpart in a hnear case. If in the equation 
(2.7) the parameter p = equation becomes hnear. If the parameter p — > 
then equation (2.7) and correspondingly equation (3.18) will be reduced to the 
hnear Black-Scholes equation but solutions (3.63)-(3.68) which we obtained 
here will be completely blown up by p — > because of the factor 1/6 = l/{up) 
in the formulas (3.63)-(3.68). This phenomena was described as well in [6], [8] 
where the solution Ur{S,t) was studied and for the complete set of invariant 
solutions of equation (3.18) with /c = 1 in [5]. 



4 Properties of solutions and peirameter-sensitivity 



We study the properties of solutions, keeping in mind that because of the 
symmetry properties (see Theorem 3.1) of the equation (2.7) we can add to 
each solution a linear function of S. 



4.1 Dependence on the constant c and terminal payoff 



First we study the dependence of the solution U2 on the arbitrary constant c. 
The constant c is the first constant of integration of the ordinary differential 
equations (3.31), (3.32). This dependence is illustrated in Figure 3. We see on 
Fig. 3 that for different values of the constant c the domains of the solutions 
are different. The dependence of the solutions M3,i, M3,2(>S', t) on the constant c 
is exemplified in Figure 4. 

We obtain a typical terminal payoff function for the solutions (3.63)-(3.68) if 

we just fix t = T. By changing the value of the constant c and by adding 
a linear function of S we arc able to modify terminal payoff function for 
the solutions. Hence we can approximate typical payoff profiles of financial 
derivatives quite well. 
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Fig. 3. Plot of solutions U2iS,t) where 
|c| = 0.01, 1., 5., 10., 20., a = 0.4, q = -4, 
b = 1.0, d = 0. The variables S lie in 
intervals S G (0., 15.) and t = 1.. 
The curves going from up to down 
with the growing value of the param- 
eter |c|. 

4-2 Dependence on time 



Fig. 4. Plot of solutions (S, t) where 
|c| = 0.01, 1., 5., 10., 20., a = 0.4, q = -4, 
b = 1.0, c? = 0. The variables S lie in 
intervals S G (0., 15.) and t = 1.. 
The curves going from up to down 
with the growing value of the param- 
eter IcI. 



All solutions depend weakly on time because of the substitution (3.49) all 
invariant solutions depends on the combination aH. As long as we take the 
volatility a to be small we obtain a dampened dependence of the solutions on 
time. 

In Fig. 5 we can see this dependence for examples of the solution ^3,1 (5*, t),U3^2{S, t). 



4.3 Dependence on the parameter p 

All solutions found in this paper have the form 

u{S,t) =w{S,t)/p, (4.69) 

where w{S,t) is a smooth function of S,t. Hence the function w{S,t) solves 
the equation (2.7) with p — 1, 



a'^S^ 'wss_ ^ . , 

+ — (l-5...,.)^ =°- 

Because of this relation any p-dependence of invariant solutions of (2.7) triv- 
ial. In particular, if the terminal conditions are fixed, u{S,T) = h{S), then 
the value u{S,t) will increase if the value of the parameter p increases. This 
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Fig. 5. Plot of solutions U3^i(5, t), U3 2(5", t) for the parameters 
a = 0.4, |c| = 10., g = —4,6 = 1.0, d = 0. The variables S lie in in- 
tervals S G (0., 15.) and t = 0.01, 1., 5., 10., 20.. 

The highest level corresponds to i = 0.01 and the lowest to i = 20.. 

dependence of hedge costs on the position of the large trader on the market 
is very natural. 



4-4 Dependence on the asset price S 



In practice one use often delta-hedging to reduce the sensitivity of a portfolio 
to the movements of an underlying asset. Hence it is important to know the 
value A defined by A = ||, where u denotes the value of the derivative 
product or of a portfolio. Using the exact formulas for the invariant solutions 
we can easily calculate A as a function of S and t. The A corresponding to 
the solution U2 is presented in Fig. 6 

The A of the solution u^^i, ^3,2 is represented on the Fig. 7. 
We see in both cases the strong dependence on S for small value of S. If 
— > 00, in both cases the A tends to a constant which is independent of time 
and the constant c. 



4-5 Asymptotic behavior of invariant solutions 



If S is large enough we have four well defined solutions Ur{S,t), Ui{S,t), 
U2{S,t), U3^2{S,t). The asymptotic behavior solutions Ur{S,t) from (3.63) and 
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Fig. 6. Plot of A for U2{S,t) and 
a = 0.4, |c| = 1., g = -4, b = 1.0, d = 0. 
The variables lie in intervals S G (0., 5.) 
and t = [0.,0.5]. 

Ui from (3.65) coincides in the main 
formula 

UriS,t), ui{S,t) r^^(3S\nS + c 



Fig. 7. Plot of the A 
for us^i{S,t),U3^2{S,t) and 

a = 0.4, |c| = l.,q = -4,b = 1.0, d = 0. 
The variables lie in intervals S G (0., 5.) 
and t G [0.01,0.5]. 

terms as S* — > oo and is given by the 
onst.S + 0{S-^/^)) (4.71) 



The exact formula for the asymptotic behavior of the function ui{S,t) for 
5* ^ cxD is given by 



u,{S, t)r.^S \og{S) + (^4 log(3) - 2 - ^ log || j - ^a' tj 

+ ^|c|e"^5-U0(5-t). (4.72) 

We see that for S ^ oo, the main term does not depend on time or on the 
value of the constant c. Moreover, this term cannot be changed by adding a 
linear function of S to the solution. 

The main terms of the solutions U2{S,t) from (3.66) and U3^2 from (3.68) 
behave similar to each other as 5 — oo; this behavior is given by the formulas 



U2{S,t) 



2^ 2\c\ 3. 
+ 36V~' 




(4.73) 



and 
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U3,2{S,t) 



s 



/2 c 3^ ^1 2^ , , 3<£i 
'-Le 32 5-4 - -_ c e 16 5 2 
3 S'^o 



o{s-i). 



(4.74) 
(4.75) 



Note that the main term in formulas (4.73) and (4.74) depends on S hn- 
early and has a coefficient which depends on the constant c only. Hence we 
can change the asymptotic behavior of the solutions U2{S,t) and ^3,1(6',^), 
1*3,2 ("S*, t) by simply adding a linear function of S to the solutions. 

For 5* in a neighborhood of 5* — > there exist just two non trivial real invari- 
ant solutions of equation (2.7), i.e. the known solution ^^(5', 1) and the new 
solution ^3,1 (5", t). Using the exact formulas for the last solution we retain the 
first term and obtain as 5" ^ for the solution 1*3,1(6', t) 



U3,i{S, t)^^ Sln{S) + {S) . (4.76) 

The main term in this formula do not depend on t or constant c and can be 
changed by adding a linear function of S to the solution. 



5 Conclusion 



In this paper we have obtained explicit solutions for the equation (2.7) and 
studied their analytical properties. These solutions are useful for a number of 
reasons. To begin with, while the payoffs of these similarity solutions cannot 
be chosen arbitrary, the payoffs can be modified using embedded constants to 
tailor a given portfolio reasonably well. For some values of the parameters c, d 
we obtain a payoff typical for futures, for other values c, d the payoff is very 
similar to the form of calls (see Figure 3 and Figure 4). Moreover, the exphcit 
solutions can be used as benchmark for different numerical methods. 

An important difference between the case of the linear Black-Scholcs equation 
and these nonlinear cases can be noticed if we consider the asymptotic for 
yS" — > 00. In the linear case the price of a Call option satisfies u{S, t) — > const- 5". 
In the nonlinear case the families of similarity solutions U2 and U3 which 
approximate the payoff of a Call option well on a finite interval [0,5], grow 
faster than linear as S* — 00; see the formulas (4.73) and (4.74) for details. 
This refiects the fact that in illiquid markets option hedging is more expensive 
than in the standard case of perfectly liquid markets. 
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